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A hysteresis model based on the assumption of fixed order magnetization reversals is pro- 
posed. The model uses one-dimensional diagram for representing states of a system despite of 
two-dimensional Preisach diagram. The distinctive feature of the model is that it is applicable to 
any system compliant with the return-point memory and includes Preisach model as a special case. 

PACS numbers: 75.60.-d, 75.60.Ej 
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INTRODUCTION 

The problem of hysteresis modeling is a subject of 
persistent interest. One of the most popular, Preisach 
model, was proposed more than 65 years ago jjj. It was 
initially based on some hypothesis concerning physical 
mechanisms of magnetization, but now it is considered 
mainly as a mathematical tool for describing various hys- 
teresis phenomena. Preisach model has many attractive 
features such as simplicity, flexibility and ability to rep- 
resent most important properties of real hysteresis sys- 
tems ||, 0|. There are known numerous applications 
of Preisach model in different areas of physics. The pos- 
tulates of Preisach model lead to the return-point mem- 
ory j^j , the property of real systems with hysteresis that 
is considered as one of the most sufficient for hysteresis 
modeling. 

In this article we intoduce a simple hysteresis model, 
which we call state vector model. It uses very different 
postulates but is similar to Preisach model in that it is 
compliant with the return-point memory and can be eas- 
ily analysed. 



DESCRIPTION OF THE MODEL 

Let us assume that states of a system can be repre- 
sented by the "state vector" x with components x a equal 
to +1 or — 1, where a = 1 ... N, and N is some large in- 
teger. 

Suppose that the magnetization M (and other physical 
values related to the system) depend on x: 



3. Each state x — {x\,X2, ■ ■ ■ , xn) of the system has 
a value of the external magnetic field 



M = M(x%, . . . , xn)- 



(1) 



The behavior of the state vector x is postulated as 
follows. 

1. When the magnetic field H increases, negative com- 
ponents x a change the sign from —1 to +1 one by 
one in order of increasing a. 

2. When the magnetic field H decreases, positive com- 
ponents x a change the sign from +1 to —1 one by 
one in order of increasing a. 



H = H(xi, . . . ,x N ) 
near which the state is stable. 



(2) 



Note that Eq. (g) must not contradict the rules 1, 2. 
While x changes according to the rules 1 and 2, the point 
defined by Eq. (|l|), (|J) form a trace on the iJM-plane. 
Due to the above postulates the system in our model 
exactly obeys the return-point memory; this will be dis- 
cussed in the next section. 

The model admits following interpretation in terms of the "en- 
semble of domains" . We may consider components indicating 
the magnetic state of the "domains" , and assume that each of them 
reverses its sign in one Barkhausen jump. In this case for the spec- 
imen of unit volume instead of Eq. (Ill) holds 



M 



rH>a X 



where m a > and x a denote the absolute value of the magnetic 
moment and its direction respectively for each domain. 

Due to the hysteresis curve symmetry, it is natural to 
assume that M{— x) = —M(x) and H(—x) — —H(x), 
which means that if the state x is stable in the external 
field H, the state —x is stable in the field — H . 

All components x a become equal to +1 in a large pos- 
itive field, and equal to —1 in a large negative field. De- 
magnetized state can be obtained as the result of ap- 
plying to the system alternating magnetic field of slowly 
decreasing magnitude, which gives the state vector with 
"disordered" components. We may take as "pure" de- 
magnetized state the pair of states (+1, — 1, +1, . . . , ±1) 
and (—1, +1, —1, . . . , ±1). The behavior of the state vec- 
tor is illustrated in Fig. [y. 

In sequel the following "continuous" formulation of the 
model will be convenient. Let us arrange all x a on the 
interval [0, 1] in the order of increasing a (Fig. ||). Func- 
tion x(£) : [0,1] — > {—1,1} determines the state of the 
system as follows: if x(£) equals to +1 (—1) on some in- 
terval, then all x a are equal to +1 (—1) on this interval. 
Instead of the rules 1 and 2, which determine the behav- 
ior of components x\ . . . , x a , we have the following rules 
for x(£). 
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FIG. 1: State vectors at different points during the magneti- 
zation process. The process starts from the demagnetized state 
(O). On the descending curve OA the magnetic field decreases, 
and x a become equal to —1 one by one from left to right. After 
the point A the magnetic field increases; starting from this point 
the components x a become equal to +1 one by one, also from left 
to right. After the point B the magnetic field decreases again and 
x a become equal to —1 in the same succession, and so on. 

1*. When the magnetic field H increases, the leftmost 
point where x(£) changes the sign from +1 to — 1, 
moves from left to right. 

2*. When the magnetic field H decreases, the leftmost 
point where x(£) changes the sign from —1 to +1, 
moves from left to right. 

To express the rules in a simple form, we have assumed 
that x(£) always changes the sign at £ = 0. With x(£) 
the right sides of Eq. ([!]) , (|2|) can be rewritten in the form 
of functionals 

M = M[x(0], H = H[x(0}. (3) 
Here H[x(£)] must not contradict 1*, 2*. 

RETURN-POINT MEMORY 

Let us show that the model can be obtained as a con- 
sequence of the return-point memory [pi. 
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FIG. 2: Illustration of the "continuous" variant of the state vector 
hysteresis model. Arrows directed up and down correspond to +1 
and —1 components of 



Return-point memory. Suppose the system is 

evolved under field H(t), where < t < T and H(0) < 
H(t) < H(T) or H(0) > H(t) > H(T). Then for a given 
initial state of the system the final state depends only on 
H(T), and is independent of the time T or the history 
H(t). 

Besides of the return-point memory the existence of 
states with some special properties is required. We sup- 
pose that a system has two states A and B, with cor- 
responding fields Ha and Hb, such that: i) the system 
evolves from the state A to the state B, if the magnetic 
field monotonically increases from Ha to Hb ; ii) the sys- 
tem evolves from the state B to the state A, if the mag- 
netic field monotonically decreases from Hb to Ha- 

In the remaining part of this section we will usually 
assume that the initial state of the system is A, all in- 
puts H(t) are continuous and piecewise monotonic and 
for each of them holds H(0) = H A , H A < H(t) < H B . 
We shall restrict our consideration to the set of states 
T,ab reachable from A by applying H (t) that satisfy the 
above conditions. 

Note that if the initial state is A, H(0) = Ha and 
H A < H(t) < H B , then H(T) = H A implies that the 
final state of the system is A, and H(T) — Hb implies 
that the final state of the system is B. 

The last statement concerning the final state B follows directly 
from the return-point memory. To ensure that it is true for the 
final state A consider the input H (t) that is applied to the system 
in the state B and is composed of H (f) before which the magnetic 
field decreases monotonically from Hg to H^. It is clear that H(t) 
and H(t) have the same final states; from the return-point memory 
follows that the final state of H (t) is A. 

Consider some input H(i), H A < H(t) < H B , H(0) = 
Ha, that is applied to the system in the state A. Let 
H(t) increase from Hq = Ha to Hi, then decrease from 
Hi to Hi, increase from Hi to H$ and so on, where 
Hi , . . . , H n —\ are successive maxima and minima oiH(t), 
and H n = H{T). 

Prime inputs. ■ Let us call a piecewise monotonic 
input H{t) prime, if H(0) = H A , H A < H(t) < H B , 
and one of the following is true: (i) the input is invari- 
able H(t) = Ha; (ii) the input H(t) is monotonically 
increasing; (iii) for all successive maxima and minima 
holds l-Hfc+i — Hk\ < \Hj~- Hk-i\ where k = 1, . . . , n — 1 
and n > 1. 
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It is reasonable do not distinguish prime inputs with 
identical values H , ■ ■ ■ , H n , because according to the 
return-point memory such inputs give the same final 
states. With this agreement any prime input is deter- 
mined by the values Hq, . . . , H n where Hq — Ha; the 
case n = corresponds to the invariable input. 

Lemma. — Each state in the set Y>ab can be obtained 
from the state A by applying a prime input. 

Proof. For any arbitrary state s in Hab exists 

piecewise monotonic input H(t) with the final state s, 
such that H(0) — Ha, Ha < H(t) < H B , < t < T. If 
H(t) = Ha at some point t', then H(t) may be replaced 
on the interval [0, t'] with invariable function H(t) = Ha- 
In the case t' = T we have the invariable input that is 
prime and has the the final state A. Otherwise, we may 
assume that H(t) does not remain constant on any subin- 
terval, because, as follows from the return-point memory, 
on each interval of monotony H(t) can be replaced with 
monotonically increasing or decreasing function without 
affecting the final state. This means that an arbitrary 
state s ^ A can be obtained by applying H(t) that 
increases from Hq = Ha to H±, then decreases from 
Hi to H2, increases from H2 to H3 and so on, where 
Ha < Hk < Hb for all k = 1 . . .n, and n is the number 
of monotony intervals. If H(t) is non-prime, then there 
are four successive extrema H^—i, Hk, Hk+i, -ff/c+2 such 
that \Hk~i — Hk\ > \Hk — Hk+i\, \H^ — Hk+i\ < \Hu+i, 
Hk+2\- According to the return-point memory, H(t) on 
the interval where H changes from H^-i to Hk+2 can 
be replaced with monotonic function without affecting 
the final state. Such replacement can be repeated not 
more than finite number of times, because the number of 
monotony intervals is finite, and finally we get a prime 
input with the final state s. The lemma is proved. 

After applying a prime input to the system in the state 
A some definite state in Y<ab is obtained, which can be as- 
signed to the prime input. According to the lemma, such 
correspondence between prime inputs and states defines 
the mapping of the set of prime inputs on the set T,ab- 

Let us consider now the set X of functions x(£) : 
[0, 1] — » { — 1, 1} which change the sign in a finite number 
of points. We can establish one-to-one correspondence 
between prime inputs and functions x(£) in the following 
way. Let us assign to the invariable input H(t) = Ha 
the function x(£) = —1. For prime inputs with n > 1 let 
us define the points 

6 = \H X - H \/\H B - H A \, 
£2 = \H 2 -H X \/\H B -H A \, 



H c 



\H n -H n _ 1 \/\H B -H A \ 



According to the definition of prime inputs we have 

<£„<... & <&<l. 

Let us assign to a given prime input the function x{^) 
that changes its sign from +1 to —1 at the point £1, 
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FIG. 3: Prime input H(t) and corresponding If H increases 

or decreases after the final point (dashed lines) x(^) changes ac- 
cording to the rules 1*, 2*. 



from —1 to +1 at the point £2 and so on. The above 
determines x(£) for any given prime input and vice versa. 
It is illustrated in Fig. || from which it is clear that x(£) 
changes according to the rules 1* and 2* of "continuous" 
state vector model, when the magnetic field decreases or 
increases after the final point of a prime input. For the 
final value of the magnetic field holds 



H = 



Ha + Hf 



Hb — Ha 



(4) 



The right side of this formula is a functional that we 
denote as H[x(£)]. As the result we have the following 
proposition. 

Proposition. The set X of functions x(£) can be 
mapped on the set of states Sab, and H [x(£)] exists, such 
that if x(£) changes according to the rules 1* and 2*, the 
state that corresponds to x(£) changes in the same way 
as the state of the system under varying external field 
H = H[x(0}. 

Once x(£) determines the state of the system, its phys- 
ical properties can be represented as the functionals with 
argument x(£); in particular, M = M[x(£)]. 

As follows from Proposition, the model is applicable to 
any system that is compliant to the return-point memory 
and has the states A and B with above mentioned prop- 
erties. The last condition is normally satisfied: we can 
take as A and B the states at the vertices of the limiting 
or any minor hysteresis loop. Note that x(£) = — 1 corre- 
sponds to the state A and x(£) = +1 corresponds to the 
state B. For the symmetric hysteresis loop Ha = —Hb 
and, according to Eq. (Q), H[x(£)] is antisymmetric; from 
the symmetry of hysteresis curves with respect to the ori- 
gin of the ifM-plane follows antisymmetry of M[x(£)]. 
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We have proved the sufficiency of the return-point 
memory for simulation by the model. The necessity also 
can be shown. Namely, it can be proved that any H [x(£)] 
which continuously and monotonically increases and de- 
creases when x(£) changes according to the rules 1* and 
2* correspondingly, is compliant with the return-point 
memory. 



FUNCTION ALS H[x(g)], M[x($)] 

Different forms of functionals H[x(£)], M[x(£)] are 
suitable for hysteresis simulation. Let us consider some 
of the simplest. 



Linear functionals 



Let us try as H[x(£)], M[x(£)] linear functionals: 



Jo 

h = [ Kt)x(t)ds, 

Jo 



(5) 
(6) 



where to(£), h(£) are some continuous positive functions. 

Consider x(£) such that = — 1 for all £ on the 

interval < £ < £i for some given £i. According to 
the model postulates it means that the magnetic field 
was decreasing before. Let the magnetic field begin to 
increase, and denote £' the point where x(£) changes the 
sign from +1 to —1, assuming £' < £1. In this case 

2 t di = \AM\, 2 / ? h(0<% = \AH\, (7) 







where AM and AH are changes of M and H starting 
from the origin of new ascending hysteresis branch; the 
same equations are true in the case of any descending 
branch. 

We can see that from Eq. (|), (§) follow that all the 
ascending hysteresis branches and all the descending hys- 
teresis branches are congruent, and can be described with 
the same function ip, which is determined by Eq. (0): 



AH = ±<p(±AM). 



(8) 



Here "+" corresponds to ascending and "— " to descend- 
ing branch. 

Let us divide the interval [0, 1] into 2n equal subinter- 
vals, and define the state x^ so that x^ ) (£) = — 1 on the 
odd subintervals and x(°\£) = +1 on the even ones. In 
accordance with the model postulates we may consider 
x (0) ( or as the demagnetized state, which is ob- 

tained via a demagnetization process performed as the 
consequence of n demagnetization cycles. From Eq. (|^), 
(||) and the continuity of m(£) and h(£) follow that 



H(x^) -> 0, M(x (0) ) -»■ when n -> oo. On the initial 
magnetization curve instead of Eq. (R) , (||) holds 

f m(0de=|Af|, f h(0d£ = \H\, 
Jo Jo 



and 



H 



±-^(±2M). 



(9) 



Two first terms of Taylor series for ip^ 1 in Eq. (JsJ) , (||) 
correspond to Rayleigh relations || ; this lead us to a con- 
clusion that the linear approximation of the functionals 
is suitable in the case of small fields. 



One nonlinear expression for H[x(£)] 

Let us consider as an example following expression 

H = a(M) f h(t)x(0<% + b(M), (10) 
Jo 

assuming that AT[x(£)] is represented by the right side of 
Eq. (||). Here a(M) and b(M) are some functions; due 
to the hysteresis loop symmetry a(M) must be even and 
b(M) must be odd. 

Prom Eq. (||), (^) we can found explicit equations for branches 
of the hysteresis curve and the initial magnetization curve 

,-Ho-&(M ) 



H = ±a(M)<p(±(M - M )) + b(M) + o(M)- 
1 



a(M ) 



H 



-a(M)tf(2M) + b(M) (M > 0). 



Here ip is determined by m(£), h(£) in the same way as it was con- 
sidered previously, and Mo , Ho denote coordinates of the beginning 
of a hysteresis branch. 

These formulae were applied for hysteresis simulation 
of low-alloyed electrical steel. Experimental and simu- 
lated curves are shown in Fig. 0. 



CONNECTION TO PREISACH MODEL 

The Preisach distribution p(h u ,h c ) describes density 
of the "domains" on the Preisach plane. Each domain 
has a shifted square hysteresis loop; h c is the coercive 
force of the domain and h u denotes the shift. Function 
p(h Ul h c ) must be positive, characterized by the syrrrrxie- 
try p(h u , h c ) = —p(—h u ,h c ), and normalized to unity. 

Let the magnetic field H always remains in the in- 
terval [— H m ,H m ], where H m is as large as desired 
but fixed; we assume also that p(h u ,h c ) is a finite 
function with support in the triangle with vertices 
(0,H m ),(0,-H m ),(H m ,0). 

Magnetization M is expressed as the integral over all 
domains 



M = M, 



dh u dh c p(h u , h c ) - M s I I dh u dh c p(h u , h c ) 
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FIG. 4: Experimental hysteresis curves of low-alloyed electrical 
steel (solid lines) and simulated curves (dashed lines). The sim- 
ulation is equally accurate in the intermediate and in the small 
fields. 



where D + and Z?_ are regions with positive and nega- 
tive domain orientation respectively. The boundary h u = 
b(h c ) between these regions is a broken line, made of seg- 
ments with positive and negative slope db(h c ) / dh c = ±1. 
The magnetization and the magnetic field are expressed 
via b(h c ): 

H m b(h c ) 

H = 6(0), M = 2M s [dh c [p(h u ,h c ) dh u . (11) 
o o 

Let us define 

x(o = —^r^bm n ), o<e<i. (12) 

H m d 4 

Function x(£) equals to ±1 and can determine the state 
of the Preisach ensemble: 



b(h c ) = H m J x(£)d£. 



(13) 




H, 



FIG. 5: A geometric illustration of the connection between 
Preisach model and the state vector model. 



Now Eq. pi] ) can be rewritten in the form of functionals 
Eq. |): 

H = H[x(Q], M = M[x(Q], (14) 



where 



and 



H[x(S)]=H m lx(Qd£, (15) 
" 



H m h c /H m 

M[x(0] = 2M S J dh c J p{h u , h c ) dh u . (16) 



The boundary b(h c ) changes in a well-known manner, 
what can be expressed in terms of This gives exactly 
the rules 1*, 2*, and in conjunction with Eq. (|lj) we 
get the state vector hysteresis model with the functionals 
of special form defined by Eq. (|l5|), (|l6|). Connections 
between two hysteresis models is illustrated in Fig. ^| 

The traditional Preisach model has some disadvantages such as 
zero valueof the turning point susceptibility and the congruency 
property n ■ Limitations of the original Preisach model take much 
weaker form in its modifications h It is possible to overcome par- 
tially the congruency property te| by introducing internal mean 
field Hmf(M), which depends on the magnetization. This exten- 
sion of Preisach model corresponds to the functionals 

H = H{x(0) - H MF (M[x(0]), M = M[x(0]- 

Another variant, the product Preisach model [^|, also can be rep- 
resented in a form of the state vector hysteresis model. In this case 
we have 

H = H[x(H)], M = G(f3B[x(£)] + M[x(0]), 

where G is "transformation function" , and the constant /3 provides 
non-zero value of the turning point susceptibility. 



h a /H„ 
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DISCUSSION AND CONCLUSIONS 

The return-point memory is the only essential con- 
dition on a system, that is necessary for simulation by 
the model. This is an obvious advantage comparative to 
Preisach model, which requires an extra condition, con- 
gruence not exhibited by majority of systems Q. 

Preisach model can be represented as a special case 
of the model, with a particular form of the function- 
als H[x(£)], M[x(£)], Explicit expressions for H[x(£)], 
M[x(£)] are found in the cases of traditional Preisach 
model and some of its extensions. From our consider- 
ation follows that a rather complex functional Eq. Jl6| ) 
that corresponds to Preisach model is not necessary for 
hysteresis simulation; it can be replaced with much sim- 
pler functionals that may not require two-dimensional 
integration (see as an example Eq. (||), ( |l0|) and Fig. ||). 

It is worth noting that a general approach for rep- 
resenting states of a system was established as a con- 
sequence of the return-point memory. This grants the 
right to use the model for any physical value y that can 
be expressed as a function of state; corresponding func- 
tional determines its behavior under varying in- 
put. Note that the model does not comprises any other 
properties except the return-point memory, and compli- 
ance with thermodynamics must be considered as a re- 



striction on the functionals. 

Different forms of the functionals could be proposed for 
different kinds of materials. Actually a class of models 
is defined, depending on a particular form of the func- 
tionals. Due to extra flexibility of the model simpler cal- 
culations and more precise simulation can be expected, 
which may assist in systematizing experimental data. 
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